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Commonly materials are classified as either electrical conductors or insulators. The theo- 
retical discovery of topological insulators (TIs) in 2005 has fundamentally challenged this 
dichotomy^. In a TI, spin-orbit interaction generates a non-trivial topology of the electronic 
band-structure dictating that its bulk is perfectly insulating, while its surface is fully con- 
ducting. The first TI candidate material put forward^ - graphene - is of limited practical 
use since its weak spin-orbit interactions produce a band-gap^ of ~0.01K. Recent reinves- 
tigation of Bi2Se3 and Bi2Te3, however, have firmly categorized these materials as strong 
three-dimensional TI's.* ^ We have synthesized the first bulk material belonging to an en- 
tirely different, weak, topological class, built from stacks of two-dimensional TI's: BiuRhsIg. 
Its Bi-Rh sheets are graphene analogs, but with a honeycomb net composed of RhBig -cubes 
rather than carbon atoms. The strong bismuth-related spin-orbit interaction renders each 
graphene-like layer a TI with a 2400K band-gap. 
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Besides being new states of matter and thus of the most fundamental scientific interest, topo- 
logical insulators (TIs) also hold promises for applications in for instance spintronics, based on the 
fact that the topological properties dictate that the metallic surface- states of TIs are spin-locked: 
theory predicts that the propagation direction of surface-electrons is robustly linked to their spin 
orientationPS. These surface- states also play a most prominent role in proposals to create Majo- 
rana fermions in microelectronic devices, the manipulation of which can be the basis for future 
topological quantum computin^^SHll, jj^ spite, of this conceptual richness, a lack of equivalent 
advances in producing new classes of TI materials has lead materials synthesis and chemistry to 
concentrate largely on further perfecting and varying upon the materials class of bismuth-based 
chalcogenides Bi2Se3 and Bi2Te3, which are confirmed three-dimensional (3D) TI^^^^ISll^ and on 
HgTe films grown with utmost care by molecular beam epitaxy which under specific conditions 
form a two-dimensional (2D) T^E. We have synthesized Bii4Rh3l9, which we will show to be the 
first member of an entirely new class of stacked 2D TIs, from a stoichiometric melt of its elements. 
The resulting thin black platelets are air- stable and can be easily cleaved. The synthetic proce- 
dure has been optimized taking into account the phase decomposition at the peritectic point of 
441°C and the strong dependence of the phase stability on the vapor-pressure. For a more detailed 
synthesis protocol we refer to the Supplementary Information (SI) . 

The crystal structure was determined by single-crystal X-ray diffraction and exhibits a pe- 
riodic alternating stacking of 2D bismuth-rhodium networks and insulating spacers, see Fig.[T]^a). 
The former, further denoted as intermetallic layer, can be understood as a decorated honeycomb 
network. It has the same hexagonal p6/mmm layer group symmetry as a graphene sheet, thus 
sharing its 24 symmetry elements. A lattice scaling of ~3.8 in Fig. [T]^b) makes this structural 
equivalence prominent. Whereas the nodes of a graphene-sheet host carbon atoms, in the inter- 
metallic layer of Bii4Rh3l9 the nodes of the net are found in the centres of triangular-prismatic 
voids (Bi-Bi: 317.5-320.6 pm). The latter appear in the hexagonal arrangement of bismuth cubes 
(Bi-Bi: 317.5-345.3 pm) which are centred by rhodium atoms (Rh-Bi: 282.0-284.9 pm). This ar- 
rangement can also be seen as a kagome-type net with the rhodium atoms of edge-sharing cubes at 
the nodes of the net. Consequently, the RhBig-cubes define the hexagonal-prismatic voids which 
are filled with iodide anions (Bi- ■ - I: 376.2-382.1 pm). As a result, an overall composition of 
(RhBi4)3l can be assigned to the intermetallic layer. Chemical-bonding analysis reveals strongly 
localized, covalent Bi-Rh bonds in the cubes and three-centred bismuth interactions in the bases of 
the triangular-prismatic voids, all together establishing a quasi-2D bimetallic network (see SI for 
details). The spacer layer consists of Bi-I zigzag chains with distorted octahedral coordination of 
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Figure 1: (a) Triclinic crystal structure of Bii4Rh3l9. Insulating layers of [Bi2l8]^^ zigzag chains 
separate the intermetallic layers [(Bi4Rh)3l]^+ which consist of hexagonal nets of edge- sharing 
RhBig cubes, (b) Honeycomb lattice of graphene scaled by a factor ~3.8 overlaid with the structure 
of the intermetallic layer. 
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bismuth(III) cations by iodide anions. Bismuth-iodine bonding between the intermetallic layer and 
the spacer (distances 376.4-455.6 pm, mean value 401.9 pm) point at weak interactions of charged 
layers and render the whole sandwiched structure of Bii4Rh3l9 = [(RhB 14)3!] ^+[31218]^^ as salt- 
hke in stacking direction. The alternate stacking of the highly-symmetric intermetallic and the 
low-symmetric spacer layer results in the reduction of the overall crystal symmetry to the triclinic 
space group PI. 

The fact that the weakly coupled intermetallic layers have the same structural symmetry 
as graphene- sheets suggests similarities in electronic structure between graphene and Bii4Rh3l9. 
Indeed a scalar-relativistic band- structure calculation for Bii4Rh3l9, where the spin-orbit coupling 
(SOC) is effectively switched off, (for details see SI), reveals the presence of two Dirac cones in 
the triclinic Brillouin-zone (BZ), which are situated at the Fermi-level, see Fig. [2^. Unfolding the 
triclinic zone to a hexagonal one, we observe that the two inequivalent Dirac cones appear very 
close to the K and K' point at the edge of the hexagonal BZ, precisely where they are in graphene 
(Fig. [2]:). The very weak hopping of electrons between the layers causes the Dirac cones to pick up 
a minor dispersion perpendicular to the plane, rendering the calculated band- structure quasi-2D. 

However, for a proper understanding of the electronic band-structure one needs go beyond 
calculations in scalar relativistic approximation. Whereas the relativistic SOC in for instance 
graphene is very weak - it is expected to open a gap on an energy scaleP of 0.01 K - bismuth 
is well-known for its strong SOC which can drive and stabilize topologically non-trivial electronic 
states. As we are thus dealing with graphene-like Dirac cones in the presence of strong SOC, the 
mechanism proposed by Kane and Mele^ to stabilize a topologically non-trivial quantum spin Hall 
state will be in action full-force in the intermetallic layers of Bii4Rh3l9. Indeed, a fully -relativistic 
band- structure calculation reveals a gapping out the Dirac cones, resulting in a calculated band-gap 
of 210 meV, corresponding to 2400 K, see Fig.[2]3. 

To confirm the topological nature of the resulting insulating state we have implemented the 
direct calculation of the four topological Z2 invariant^^IZHIll i/^; (i/^^ 1^2, ^'3) in the FPLO band- 
structure codeP' (see SI for details). This calculation is based on an analysis of the wavefunction 
parity eigenvalues at the eight time-reversal invariant points of the band-structure, as illustrated in 
Fig. |2]F. Because of the stacking of the quasi-2D intermetallic planes in which SOC has gapped 
out the Dirac cones, one expects Bii4Rh3l9 to be a weak TI and indeed we find uq = 0. The other 
topological invariants are calculated to be {ui, 1/2, z/3) = (0, 0, 1). The fact that uq = and = I 
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without spin-orbit interaction 



with spin-orbit interaction 




Figure 2: (a) Scalar relativistic band-structure of triclinic Bii4Rh3l9. On the line FX a Dirac cone 
is present, (b) Fully relativistic band-structure where spin-orbit interaction opens up a gap of ~210 
meV (shaded area), (c) Positions of Dirac points (DP) in the basal-plane of the triclinic BZ (dashed 
hnes) and the hexagonal BZ (full lines). DPs are on the line FK in the triclinic BZ, corresponding 
to a position close to the K points of the unfolded hexagonal BZ. (d) Illustration of the Dirac cones 
in the situation without spin-orbit interaction and (e) the topological gap induced by spin-orbit 
interaction, (f) Parity eigenvalues at the 8 time reversal invariant points in the 3D BZ, where 
corresponds to the direction perpendicular to the Bi-Rh planes. The parity eigenvalues lead to the 
topological invariants z/q; (i^i, ^2, 1^3) = 0; (0, 0, 1). Projected parities on the planes perpendicular 
to gi and g2 illustrate the presence of topological edge-states on these surfaces. 
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proves that Bii4Rh3l9 is a weak TI - the first synthesized material in this topological class - and 
confirms that the intermetallic planes form sheets of Quantum Spin Hall states that are stacked 
along the c-axis. It implies that in the Altland-Zimbauer classification^' Bii4Rh3l9 belongs to the 
symplectic (All) topology class. On surfaces perpendicular to (0,0, 1) - and therefore parallel 
to the intermetallic planes, which are the natural cleaving planes of the material - topological 
surface-states will thus be absent. At any other surface an even number of Dirac cones appears, 
each having a strongly anisotropic group velocity due to the quasi-2D nature of the bulk band- 
structure. The metallicity of these surface-states is stable against disorder, which, as for strong TIs 
in the symplectic class, does not act as a source of localizatiorP^Hl 

The band-structure obtained from the electronic structure calculations can be tested experi- 
mentally by Angle Resolved Photoemission (ARPES), which in particular can provide the experi- 
mental value of the electronic band-gap, if by slight electron doping the Fermi-level of the material 
is pushed into the conduction band. The ARPES spectra in Fig. |3^,d, show that indeed the ma- 
terial is n-doped, which is consistent with a slight iodine-deficiency related to its vapor-pressure 
or alternative causes for charge redistribution at the surface probed by ARPES. It offers a unique 
opportunity to check the requisites of the topological scenario experimentally, in a full analogy 
with for instance the well-known caseP^of Bi2Se3. 

We compare the ARPES results with the fully relativistic band-structure calculations in the 
unfolded BZ selecting the FK high-symmetry direction (Fig. [3^,b, for other cuts see SI). This is 
a direction where the Dirac cones are found in scalar-relativistic approximation. The agreement 
between experimental data and fuUy-relativistic calculations is remarkable. First of all, the ob- 
served gap, as determined considering the distance between the features of the integrated (over 
three BZ) spectral weight, is ~270 meV, which is consistent with the relativistic band-structure. 
Also the number and behavior of the dispersing features clearly seen in ARPES intensity plots 
are essentially captured by the bands projected to the hexagonal BZ. The prominent examples are 
the inequivalence of F and r2 points and structures between K and K2, both also present in the 
calculated band-structure. 

This agreement confirms experimentally the correctness of the calculated band-structure in 
presence of spin-orbit coupling - the bands for which we have calculated the topological invariants 
to be uq; (z^i, z/2, z/3) = 0; (0,0,1). This topologically non-trivial state we have understood in 
terms of the graphene-like structure of the intermetallic Bi-Rh planes, for which two Dirac cones 
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Figure 3: Band- structure as measured by angle-resolved photo-emissiorP^ compared to the calcu- 
lated one. A small amount of electron doping brings the Fermi-level into the conduction band, so 
that the topological band-gap becomes clearly visible: (a) Momentum-energy intensity plot along 
the high-symmetry FK direction of the unfolded Brillouin zone, (b) Unfolded fully relativistic 
band-structure along the same direction, (c) Momentum distribution at 0.4 eV binding energy in- 
tegrated within 50 meV window. White dashed lines show the trigonal, surface projected BZ. (d) 
Typical momentum-energy distribution taken along the cut which corresponds to = in panel 
(c). 
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are present at the Fermi-level in a calculation without spin-orbit coupling. Experimentally the 
spin-orbit coupling is of course unavoidable and switching it on in the calculation gaps out the 
Dirac cones and generates the weak topological 0; (0, 0, 1) state. As a consequence spin-locked 
topological surface states will be present at any crystal face that is not parallel to the (001) plane. 
Such faces, however, do not correspond to natural cuts of the crystal and the associated surface 
roughness has prevented us so far from observing these spin-polarized surface states by ARPES 
on, for instance, (100) surfaces. In principle (spin-polarized) scanning tunneling microscopy and 
spectroscopy should be able to directly probe these states at step-edges of the natural (001) cleaving 
plane. 

Finally, the momentum distribution of ARPES intensity at 400 meV binding energy, just 
below the calculated Fermi level of Bii4Rh3l9, clearly shows a hexagonal pattern (Fig. [3]:), as is 
expected from the band- structure calculations where the top of the valence band is formed by a 
dispersive feature that is relatively flat between the K points. The hexagonal pattern justifies the 
unfolding of the BZ carried out to facilitate the comparison with ARPES. At the same time its 
slight irregularity reflects the over-all triclinic crystal symmetry. 

The hexagons observed in ARPES emphasize the structural and electronic similarities of 
graphene and Bii4Rh3l9, which again it shares with the chemically and structurally closely re- 
lated compoun(P3^ BiiaPtsIy, for instance. The compelling difference with graphene is that the 
large spin-orbit coupling drives Bii4Rh3l9 electronically into a topologically insulating state, cor- 
responding to a 3D stack of 2D quantum spin Hall states, which is very different from the strong 
TI states observed so far in bulk materials such as Bi2Se3/Bi2Te3 and which is predicted to leave 
marked signatures on electronic transport through its spin-polarized, stacked quasi- ID topological 
surface-state j^^ l ^^ l. 
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